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Abstract 

Using the methods developed for the Bianchi I case we have shown that a boostrap argu- 
ment is also suitable to treat the future non-linear stability for reflection symmetric solutions 
of the Einstein-Vlasov system of Bianchi types II and VIo. These solutions are asymptotic 
to the Collins- Stewart solution with dust and the Ellis-MacCallum solution respectively. We 
have thus generalized the results obtained by Rendall and Uggla in the case of locally rota- 
tionally symmetric Bianchi II spacetimes to the reflection symmetric case. However we needed 
to assume small data. For Bianchi VIo there is no analogous previous result. 

Contents 

1 Introduction 

2 Relativistic kinetic theory 

2.1 The Einstein-Vlasov system 

2.2 Energy momentum tensor and characteristics 

2.3 The initial value problem 

3 Bianchi spacetimes 

3.1 Definition of Bianchi spacetimes 

3.2 Description of Bianchi spacetimes via the metric approach 

3.3 3+1 Decomposition of the Einstein equations 

3.4 Time origin choice and new variables 

3.5 Vlasov equation with Bianchi symmetry 

4 Bianchi A spacetimes 

4.1 Definition and classification of Bianchi A spacetimes 

4.2 Reflection symmetry and vanishing tilt 

4.3 Some formulas for the diagonal case 

4.3.1 Expressions for diagonal Bianchi II 

4.3.2 Expressions for diagonal Bianchi VIo 

5 Special solutions 

5.1 The Kasner solution and generalized Kasner exponents 

5.2 The Collins- Stewart solution 

5.3 The Ellis-MacCallum solution 



i 



Einstein-dust with small data 

6.1 Einstein-dust system 

6.2 Linearization of Einstein-dust around Collins-Stewart . 

6.3 Linearization of Einstein-dust around Ellis-MacCallum 

The bootstrap argument 

7.1 Bootstrap assumptions 

7.1.1 Bootstrap assumptions for Bianchi II 

7.1.2 Bootstrap assumptions for Bianchi VIo 

7.2 Estimate of the mean curvature 

7.3 Estimate of the metric 

7.4 Estimate of P 

7.5 Closing Bianchi II 

7.5.1 Estimate for £_ 

7.5.2 Bootstrap assumptions for the other time variable 

7.5.3 Estimate for E + and N\ 

7.6 Closing Bianchi VIo 

7.7 Results of the bootstrap argument 

Main results 

8.1 Arzela Ascoli 

8.2 Optimal estimates 

8.3 Kasner exponents 

8.4 Estimates of the energy momentum tensor 

Conclusions and Outlook 



ii 



1 Introduction 



A starting point to understand general cosmological models are the homogeneous models. There 
are a lot of results concerning this subject and in particular two books [34], [39] which are an 
excellent introduction and a great summary of many of the results obtained. 

In general the focus has been on the fluid model since it appears (theoretically) relatively 
natural when dealing with isotropic universes and from observations we also know that the Universe 
is almost isotropic. However to have a deeper understanding of the dynamics one should go beyond 
the study of isotropic universes. General statements may vary then depending on the choice of 
the matter model. It is also important to note as is pointed out in [pj that a quasi-isotropic epoch 
is compatible with all Bianchi models and thus it is interesting to study the dynamics of all the 
different types. 

We will deal with the future asymptotics of some homogeneous cosmological models within the 
so called Bianchi class A and the matter is described via an ensemble of free falling particles also 
called collisionless matter. For all the models treated here the fundamental questions are on a 
firm ground, i.e. future geodesic completeness has been shown for these models [53], |24| . 

Concerning the late time behaviour of the Universe one believes in the cosmic no hair con- 
jecture. This conjecture states roughly speaking that all expanding cosmological models with a 
positive cosmological constant approach asymptotically the de Sitter solution. For the Einstein- 
Vlasov system isotropization could be shown for all non-type IX Bianchi cosmologies |19| . 

In absence of a cosmological constant there are also different results concerning the future. 
For the Einstein- non- linear scalar field system we refer to [31], [35], pQ and [7]. The late-time 
behaviour of Bianchi spacetimes with a non-tilted fluid is well understood [3H] , [H] . In particular 
all non-tilted perfect fluid orthogonal Bianchi models except IX with a linear equation of state 
where < 7 < | are future asymptotic to the flat FL model [TB] . For other values of 7 one cannot 
expect isotropization for most of the Bianchi models. However there is an important characteristic 
of the future asymptotics for the Bianchi types I, II and VIo, namely that the spacetimes considered 
tend to special (self-similar) solutions. For expanding models it is reasonable to expect that the 
dispersion of the velocities of the particles will decay. The conjecture now is that for the Einstein- 
Vlasov system the spacetimes of certain Bianchi type will tend to the same special solution as for 
the corresponding Einstein-dust system. 

This has been achieved already for locally rotationally symmetric (LRS) models in the cases 
of Bianchi I, II, III, VIII and IX ^ , [3D] , [55] , Q3], [3J and for reflection symmetric models in the 
case of Bianchi I [5lj]. These results have been obtained using dynamical systems theory. 

The Einstein- Vlasov system remains a system of partial differential equations (PDE's) even if 
one assumes spatial homogeneity. The reason is that although the distribution function written 
in a suitable frame will not depend on the spatial point, the dependence with respect to the 
momenta remains (since the Vlasov equation is defined on the mass shell) . However in the results 
mentioned a reduction to a system of ordinary differential equations was possible due to the 
additional symmetry assumptions. This is no longer possible if one drops some of these additional 
symmetries (see [50] for the reasons). Thus if one wants to generalize these results the theory of 
finite dimensional dynamical systems is not enough. 

Most of the results obtained until now rely on the theory of dynamical systems. Thus one 
might be tempted to use techniques coming from the theory of infinite-dimensional dynamical 
systems. The first important difficulty would be to choose the suitable (weighted) norm. Another 
one is that important theorems which have been used for the finite-dimensional case cannot be 
used here. All this may work, but this is not the approach taken here. 

Here the main tool used is a bootstrap argument which is often used in non-linear PDE's. We 
will present results concerning the late-time behaviour of some expanding Bianchi A spacetimes 
with collisionless matter where we have assumed small data. This assumption will be specified 
later, but roughly it means that the universe is close to the special self-similar solution mentioned 
earlier and that the velocity dispersion of the particles is small. 

The results obtained are as follows. For reflection symmetric Bianchi II and reflection sym- 
metric Bianchi VIq we have been able to show that their late-time behaviour remains the same if 
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the LRS condition is dropped. We will show that these spacetimes, reflection symmetric Bianchi 
II and reflection symmetric Bianchi VIo, will tend to solutions which are even more symmetric. 
In the case of Bianchi II we will show that it will become LRS, a Bianchi model whose isometry 
group of the spatial metric is four-dimensional. In this case there exists a one-dimensional isotropy 
group and one can show that a spacetime of Bianchi class A admits a four-dimensional isometry 
group, if and only if two structure constants are equal and if the corresponding metric components 
are equal as well. Bianchi VIo cannot be LRS, however it is compatible with an additional discrete 
symmetry (Appendix B.l of [4]). The analysis of the asymptotics shows that the Bianchi VIo 
spacetimes tend to this special class. Note that for VIo there is no corresponding LRS/previous 
result. 

To be more precise we show that the reflection symmetric solutions of Bianchi type II and 
VIo are asymptotic to the Collins-Stewart solution with dust and the Ellis-MacCallum solution 
respectively. The asymptotic behaviour of the metric and the energy-momentum tensor is studied 
in detail. 

All the results show that the dust model usually assumed in observational cosmology in the 
'matter-dominated' Era is robust. Another way of saying the same is that asymptotically colli- 
sionless matter is well approximated by the dust system. 

The paper is organized as follows. In the following section we will present the general basic 
equations, namely the Einstein- Vlasov system. In sections three and four we explain the symmetry 
assumptions and deduce the corresponding equations of the Einstein- Vlasov system which are then 
summarized in section five. Afterwards we present some special solutions and their linear stability 
as a pre-stage of our main argument in section eight: the bootstrap argument. This argument 
is refined in section nine and leads to our main results. In our last section we have a discussion 
about our results and present some possible future directions. 

2 Relativistic kinetic theory 
2.1 The Einstein- Vlasov system 

A cosmological model represents a universe at a certain averaging scale. It is described via a 
Lorentzian metric g a p (we will use signature - + + +) on a manifold M and a family of fundamental 
observers. The metric is assumed to be time-orientable, which means that at each point of M the 
two halves of the light cone can be labelled past and future in a way which varies continuously 
from point to point. This enables to distinguish between future-pointing and past-pointing timelike 
vectors. This is a physically reasonable assumption from both a macroscopic point of view e.g. the 
increase of entropy and also from a microscopic point of view e.g. the kaon decay. As we already 
mentioned in the introduction one has also to specify the matter model and this we will do in this 
section. The interaction between the geometry and the matter is described by the Einstein field 
equations (we use geometrized units, i.e. the gravitational constant G and the speed of light in 
vacuum c are set equal to one): 

G a f) = 87rT Q/ 3 

where G a p is the Einstein tensor and T a p is the energy-momentum tensor. For the matter model 
we will take the point of view of kinetic theory [35]. The sign conventions of [5S| are used. Also 
the Einstein summation convention that repeated indices are to be summed over is used. Latin 
indices run from one to three and Greek ones from zero to three. 

Consider a particle with non-zero rest mass which moves under the influence of the gravitational 
field. The mean field will be described now by the metric and the components of the metric 
connection. The wordline x a of a particle is a timelike curve in spacetime. The unit future- 
pointing tangent vector to this curve is the 4- velocity v a and p a = mv a is the 4-momentum of the 
particle. Let T x be the tangent space at a point x a in the spacetime M, then we define the phase 
space for particles of arbitrary rest masses P to be the following set: 

P = {(x a ,p a ) : x a e M, p a G T x , p aP a < 0, p° > 0} 
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which is a subset of the tangent bundle TM = {(x a ,p a ) : x a £ M, p a £ T x }. For particles of the 
same type and with the same rest mass to which is given by the mass shell relation: 

PaP a = -m 2 

we have the phase space P m for particles of mass to: 

P m = {(x a ,p a ) : / e M, p a £ T x , PaP a = -to 2 , p° > 0} 

We will consider from now on that all the particles have equal mass to. For how this relates to the 
general case of different masses see [3]. We will choose units such that m = 1 which means that 
a distinction between velocities and momenta is not necessary. We have then that the possible 
values for the 4-momenta are all future pointing unit timelike vectors. These values form the 
hypersurface: 

Pi. = {(x a ,p a ) : x a £ M, p a £ T x , PaP a = -1, p a > 0} 

which we will call the mass shell. The collection of particles (galaxies or clusters of galaxies) 
will be described (statistically) by a non- negative real valued distribution function f(x a ,p a ) on 
Pi . This function represents the density of particles at a given spacetime point with given four- 
momentum. A free particle travels along a geodesic. Consider now a future-directed timelike 
geodesic parametrized by proper time s. The tangent vector is then at any time future-pointing 
unit timelike. Thus the geodesic has a natural lift to a curve on Pi by taking its position and 
tangent vector. The equations of motion thus define a flow on Pi which is generated by a vector 
field L which is called geodesic spray or Liouville operator. The geodesic equations are: 



dx a dp a ra « 7 

-ds~ =P 5 -ds- = - T ^ P P 



where the components of the metric connection, i.e. T a p^ = g(e a , V 7 ep) = g a s^fj 7 can be ex- 
pressed in the vector basis e a as [(1.10.3) of |36|]: 



r Q/ 3 7 = ^[epida-y) + e i(9i3a) + e a {g lfi ) + rfpgaS + V 5 ai 9i3S - Vp a 9is] (1) 
The commutator of the vectors e a can be expressed with the following formula: 

r i 7 

l e oi e P\ — Va^l 

where ry a o are called commutation functions. 

The restriction of the Liouville operator to the mass shell is defined as: 

dx a d dp a d 
— _i_ 

ds dx a ds dp a 
Using the geodesic equations it has the following form 

L - p dx* P dp- 



This operator is sometimes also called geodesic spray. If we denote now the phase space density 
of collisions by C(/), then the Boltzmann equation in curved spacetime in our notation looks as 
follows: 

L(f) = C(f) 

describes the evolution of the distribution function. Between collisions the particles follow geodesies. 
We will consider the collisionless case which is described via the Vlasov equation: 

L(f) = 
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2.2 Energy momentum tensor and characteristics 



The unknowns of our system are a 4-manifold M, a Lorentz metric g a p on this manifold and the 
distribution function / on the mass shell Pi defined by the metric. We have the Vlasov equation 
defined by the metric for the distribution function and the Einstein equations. It remains to define 
the energy-momentum tensor T a p in terms of the distribution and the metric. Before that we need 
a Lorentz invariant volume element on the mass shell. A point of a the tangent space has the 
volume element \g^\ ? dp dp 1 dp 2 dp 3 {g^ is the determinant of the spacetime metric) which is 
Lorentz invariant. Now considering p° as a dependent variable the induced (Riemannian) volume 
of the mass shell considered as a hypersurface in the tangent space at that point is 

w = 2H(p a )6(p aP a + m 2 )\g^\Up°dp 1 dp 2 dp 3 

where S is the Dirac distribution function and H{p a ) is defined to be one if p a is future directed 
and zero otherwise. We can write this explicitly as: 

™ = \p \- 1 \g {4) \3dp 1 dp 2 dp 3 

Now we define the energy momentum tensor as follows: 

T a p = J f{x a 1 p a )p a p p zu 

One can show that T a p is divergence-free and thus it is compatible with the Einstein equations. 
For collisionless matter all the energy conditions hold (for details we refer to [27] ). In particular 
the dominant energy condition is equivalent to the statement that in any orthonormal basis the 
energy density dominates the other components of T a p, i.e. T Q( g < Too for each a, (3 (P. 91 of |12|V 
Using the mass shell relation one can see that this holds for collisionless matter. The non-negative 
sum pressures condition is in our case equivalent to g a bT ab > 0. 

The Vlasov equation in a fixed spacetime can be solved by the method of characteristics: 

dXa _ pa. dpa _ r a p p pl 

Let X a (s,x a ,p a ), P a (s,x a ,p a ) be the unique solution of that equation with initial conditions 
X a (t, x a ,p a ) = x a and P a (t,x a ,p a ) = p a . Then the solution of the Vlasov equation can be 
written as: 



f{x a ,p a ) = f (X a (0,x a ,p a ),P a (0,x a ,p a )) 

where fo is the restriction of / to the hypersurface t = 0. It follows that if fo is bounded the same 
is true for /. We will assume that / has compact support in momentum space for each fixed t. 
This property holds if the initial datum fo has compact support and if each hypersurface t = to 
is a Cauchy hypersurface. Note that there is no obvious reason why a solution of the Boltzmann 
equation with compactly supported initial data should have compact support |27) . 



2.3 The initial value problem 

Before coming to our symmetry assumption we want to briefly introduce the initial value problem 
for the Einstein- Vlasov system. For a general introduction to the initial value problem in general 
relativity we refer to J3Hj and for the Einstein- Vlasov system in particular we refer to [37] . In 
general the initial data for the Einstein-matter equations consist of a metric g a b on the initial 
hypersurface, the second fundamental form k a b on that hypersurface and some matter data. Thus 
we have a Riemannian metric g a b, a symmetric tensor k a b and some matter fields defined on an 
abstract 3-dimensional manifold S. 

Solving the initial value problem means embedding S into a 4-dimensional M on which are 
defined a Lorentzian metric g a p and matter fields such that g a b and k a b are the pullbacks to S of 
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the induced metric and second fundamental form of the image of the embedding of S while / is 
the pullback of the matter fields. Finally g a @ and / have to satisfy the Einstein-matter equations. 

For the Einstein- Vlasov system it has been shown [S] that given an initial data set there exists 
a corresponding solution of the Einstein- Vlasov system and that this solution is locally unique up 
to diffeomorphism (see also theorem 1.1 of [37]). The extension to a global theorem has not been 
achieved yet. However if one assumes that the initial data have certain symmetry, this symmetry 
is inherited by the corresponding solutions (see 5.6 of |10| for a discussion). In particular for the 
case we will deal with, i.e. expanding Bianchi models (except type IX) coupled to dust or to 
collisionless matter the spacetime is future complete (theorem 2.1 of 1241). 

3 Bianchi spacetimes 

3.1 Definition of Bianchi spacetimes 

We start with the definition of homogeneity of spacetimes taken from chapter 5.1 of [40]. 

Definition 1. A spacetime (M,g a p) is said to be (spatially) homogeneous if there exists a one- 
parameter family of spacelike hypersurfaces St foliating the spacetime such that for each t and for 
any points P,Q € St there exists an isometry of the spacetime metric, g a p, which takes P into Q. 

The basis for the classification of homogeneous spacetimes is the work of Bianchi [2] which was 
introduced to cosmology by Taub |37j . Here we will use the modern terminology and we define 
Bianchi spacetimes as follows: 

Definition 2. A Bianchi spacetime is defined to be a spatially homogeneous spacetime whose isom- 
etry group possesses a three-dimensional subgroup G that acts simply transitively on the spacelike 
orbits. 

Not all homogeneous spacetimes are Bianchi spacetimes. But the only case where G does not 
act simply transitively or does not possess a subgroup with simply transitive action are the so 
called Kantowski-Sachs models. The Bianchi models can be subclassified into two classes [8]: class 
A and B. Later we will only deal with Bianchi class A, however all the equations in this s are valid 
for Bianchi spacetimes in general. 

The only Bianchi spacetimes which admit a compact Cauchy hypersurface are Bianchi I and 
IX. In order to be not that restrictive we will consider locally spatially homogeneous spacetimes. 
They are defined as follows. Consider an initial data set on a three-dimensional manifold 5*. Then 
this initial data set is called locally spatially homogeneous if the naturally associated data set on 
the universal covering S is homogeneous. For Bianchi models the universal covering space S can 
be identified with its Lie group G (see [23], [2J] for details). 

3.2 Description of Bianchi spacetimes via the metric approach 

A Bianchi spacetime admits a Lie algebra of Killing vector fields with basis ki , k2 , k3 and structure 
constants C^ b , such that: 

[k a ,k b ] = -c c ab k c . 

The Killing vector fields k a are tangent to the group orbits which are called surfaces of homogeneity. 
If one chooses a unit vector field n normal to the group orbits we have a natural choice for the time 
coordinate t such that the group orbits are given by a constant t. This unit normal is invariant 
under the group, i.e: 

[n, k a ] = 

One can now choose a triad of spacelike vectors e a that are tangent to the group orbits: 

g(n,e a ) = 
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and that commute with the Killing vector fields: 



[e a ,k b ] = 

A frame {n, e a } chosen in this way is called a left invariant frame and it is generated by the right 
invariant Killing vector fields. Since n is hypersurface orthogonal the vector fields e a generate a 
Lie algebra with structure constants if ah . It can be shown that this Lie algebra is in fact equivalent 
to the Lie algebra of the Killing vector fields. Thus one can classify the Bianchi spacetimes using 
either the structure constants or the spatial commutation functions of the basis vectors. The 
remaining freedom in the choice of the frame is a time-dependent linear transformation, which can 
be used to introduce a set of time- independent spatial vectors E a : 

[E a ,n] = 0. 

The corresponding commutation functions are then constant in time and one can make them equal 
to the structure constants: 

[E Q ,E b ] = Cq 6 E c 

This is our choice which is sometimes called the metric approach. If W a denote the 1-forms dual 
to the frame vectors E a the metric of a Bianchi spacetime takes the form: 

i g=-dt 2 + g ab (t)W a W b (2) 

where g a b (and all other tensors) on G will be described in terms of the frame components of the 
left invariant frame which has been introduced. A dot above a letter will denote a derivative with 
respect to the cosmological time t. 

3.3 3+1 Decomposition of the Einstein equations 

We will use the 3+1 decomposition of the Einstein equations as made in j28J. Comparing our 
metric with (2.28) of [55] we have that a = 1 and (3 a = which means that the lapse function 
is the identity and the shift vector vanishes. There the abstract index notation is used. We can 
interpret the quantities as being frame components. For details we refer to chapter 2.3 of |28| . 
There are different projections of the energy momentum tensor which are important 

p = T 00 

J a — -La 
S a b = T a b 

where p is the energy density and j a is the matter current. 
The second fundamental form k a b can be expressed as: 

9ab = -2k a b- (3) 

The Einstein equations: 

Kb = Rab + k kab ~ ^k ac kl - 8n(S a b - \gabS) - ^PQab (4) 

where we have used the notations S = g ab S a b, k = g ab k a b, and R a b is the Ricci tensor of the three- 
dimensional metric. The evolution equation for the mixed version of the second fundamental form 
is (2.35) of EFJ: 

kt = R a b +k kt - 87rS£ + 4tt6£ {S - p) (5) 
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From the constraint equations since k only depends on the time variable we have that: 



R-k ab k ab + k 2 = 16717? (6) 
V a fc ab = 87rj b (7) 

where R is the Ricci scalar curvature. 

Another useful relation concerns the determinant g of the induced metric ((2.30) of |28 ): 

j t {\ogg) = -2k (8) 

Taking the trace of (J5J) : 

k = R + k 2 + 4ttS - 12np (9) 
With ^ one can eliminate the energy density and © reads: 

k= ^(k 2 + R + 3k ab k ab ) + 4irS (10) 

Finally if one substitutes for the Ricci scalar with ([5]): 

k = k ab k ab + 4tt(5 + p) (11) 



3.4 Time origin choice and new variables 

Now with the 3+1 formulation our initial data are (gij (to), kij (to), f(to)), i.e. a Riemannian metric, 
a second fundamental form and the distribution function of the Vlasov equation, respectively, on 
a three-dimensional manifold S(to). This is the initial data set at t = to for the Einstein- Vlasov 
system. 

We assume that k < for all time following [S3] (see comments below lemma 2.2 of [23J). 
This enables us given k(to) to set without loss of generality to = —2/k(to), since future geodesic 
completeness is known and one can always add an arbitrary constant to the time origin. The 
reason for this choice will become clear later and is of technical nature. 

We will now introduce several new variables in order to use the ones which are common in 
Bianchi cosmologies and to be able to compare results. We can decompose the second fundamental 
form introducing a ab as the trace-free part: 

k a b = o ab - Hg ab (12) 



k ab k ab - o ab o ab + 3H 2 (13) 



Using the Hubble parameter: 



we define: 



and 



H = —k 
3 



K = f (14) 



S + = -i(S^+Sl) (15) 

e_ = -^=(e2-e£) (i6) 
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Define also: 

Q = 8np/3H 2 (17) 
*=- l -H2 ( 18 ) 

The time variable r is dimensionless and sometimes very useful. From (O we obtain the constraint 
equation: 

J^(R-a ab a°») = n-l 
and from (|10p the evolution equation for the Hubble variable: 



Combining the last two equations with ([5J we obtain the evolution equations for £_ and E 



+ • 

£+ = g[ 2fl " 3 ^ + Jg) - S + (3 + f,) + ^(35 2 2 + 351 - 25)] (21) 

S - - ^ [ " ( + ^ + ] (22) 

3.5 Vlasov equation with Bianchi symmetry 

Since we use a left-invariant frame / will not depend on x a and the Vlasov equation takes the 
form: 

P dt P dp- ~° 

It turns out that the equation simplifies if we express / in terms of Pi instead of p l what we can 
do due to the mass shell relation: 

dt dp a 

Because of our special choice of frame the metric has the simple form ^j. This has the consequence 
that only the spatial components of the metric connection remain and that the first three terms 
of ([1} vanish. Due to the fact that we are contracting and the antisymmetry of the structure 
constant we finally arrive at: 

ff + W-'Ctlft-g-O (23) 
From (|2"3")l it is also possible to define the characteristic curve V a : 

{V^CIV'V, (24) 



dt 

for each Vi (t) = vi given t. Note that if we define: 

V = g ij ViV 3 (25 
due to the antisymmetry of the structure constants we have with (|24D : 



w = £<""> v ' v ' (26) 
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Let us also write down the components of the energy momentum tensor in our frame: 



Too = / /(^VvW^P 3 (27) 
T j = -J f(t, p a )p 1 ^-gdp 1 dp 2 dp i (28) 
Tij = [ f{t 1 p a )p l p 3 {p*)- 1 J-gdp 1 dp 2 dp i (29) 



4 Bianchi A spacetimes 

4.1 Definition and classification of Bianchi A spacetimes 

In the last chapter we have presented the Einstein- Vlasov system with Bianchi symmetry. However 
our results concern only a special class of the Bianchi spacetimes, namely that of class A. 

Definition 3. A Bianchi A spacetime is a Bianchi spacetime whose three-dimensional Lie algebra 
has traceless structure constants, i.e. C% = 0. 

In that case there is a unique symmetric matrix, called commutator matrix with components 
v % i such that the structure constants can be written as follows (lemma 19.3 of |33|): 

C% c = e hcd v da (30) 

The transformation rule of the commutator matrix under a change of basis of the Lie algebra 
can be used to classify the Bianchi class A Lie algebras. It is possible to diagonalize v and the 
diagonal elements v\ , v 2 and 1/3 can be used to classify the different Bianchi types of class A (for 
details see chapter 19.1 of )• We will study Bianchi II and VIo- For Bianchi II we have that 
the only non- vanishing element of v is v\ = 1, thus from (|30j) we see that for Bianchi II the only 
non-vanishing structure constants are: 

C\ z = 1 = -C\ 2 (31) 

In the case of Bianchi VIo the only non- vanishing elements of v are 1/2 = 1 and v-$ = — 1 thus the 
only non- vanishing structure constants are: 

C 31 = 1 = — C 13 , C21 = 1 = — C 3 2 (32) 



4.2 Reflection symmetry and vanishing tilt 

We will assume an additional symmetry namely the reflection symmetry such that the matter 
current vanishes. This will be a restriction to the diagonal case. The reflection symmetry has 
been defined in (2.10) of [25] for the case of Bianchi I, but one can define this for other Bianchi 
types as well with the difference that the distribution function now will depend in general on the 
time variable: 

f(t,Pl,P2,P3) = f(t,-pi,~P2,P3) = f{t,Pi,-p2,~P3) 

Suppose that the distribution function is initially reflection symmetric and the metric and the 
second fundamental form are initially diagonal. Then from (|28p - (|2"9"|) we see that the energy- 
momentum tensor is diagonal as well. From @ and ^ we can see that the metric and the second 
fundamental will remain diagonal. This symmetry implies in particular that there is no matter 
current, which means that there is no 'tilt'. Initial data with symmetry lead to solutions of the 
Einstein- Vlasov equations with symmetry, for a proof of this fact we refer to 9.2-9.3 of |28) . 
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4.3 Some formulas for the diagonal case 

Now we will introduce formulas which are valid in the diagonal case. There (ijk) denotes a 
cyclic permutation of (123) and the Einstein summation convention is suspended for the first two 
formulas. Let us define: 



9i; 

m = Vi 



9jj9kk 



The Ricci tensor is given by (11a) of [1]: 



We define: 



In the diagonal case we have: 



H 



1 O-afcCT 



6 H 2 

from which follows that the constraint equation can be written in the following way: 

Y? + + £ 2 _ = n - 1 - -j^ R (33) 
Now we proceed to use them for the cases of Bianchi II and VIo . 

4.3.1 Expressions for diagonal Bianchi II 

For Bianchi II we have then: 

R\ = -Rl = -R\ = -R= \n\ 



From the constraint equation (|33p we obtain: 

+ Y?_ = 1 - n - ^-Nl 
and from ([2U)) we obtain the equation for the evolution of H: 



„ , , N 3 iV 2 3 , 9 , s 4nS , s 

) = 2-2l + 2^ +E -) + 3^ (34) 



and the evolution equation for S + and £_: 

S+ = H[±N* - (3 + J^)£+ + ^(S 2 + 5f - 25?)] (35) 
S- = H[-(3 + ^)£_ + ^(^ 2 - 5 3 3 )] (36) 

From the definition @ for the second fundamental form the evolution equation for n\ follows: 

d 



In terms of Nf: 



j- t {n\) = 2(-Aa + -H)n\ 



j t {Nl) = -2NlH{^ + + l + ±L) 
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or 



Ni = -iVifr(4E+ + 1 + -p) (37) 
Let us write the equations (|35l) - (|37l) with r and g: 

S+ = l^f - (2 - g )E+ + ^ (5| + Si - 2Sl) (38) 

4-7T 

EL = -(2 - g)E_ + ~^{Sl Si) (39) 
iV( =N 1 (q-4S+) (40) 
Note that these equations are the same as (6.21) of [35] with 7 = 1 if one sets S = in 

4.3.2 Expressions for diagonal Bianchi VIo 

For Bianchi VIo we have then: 

R\=R=-l{n 2 -nz) 2 



fla = -fli = o («£-»§) 



The constraint equation ([55)) is: 



e^ + eL = i-fi-^(iv 2 -iv 3 ) 2 



and the evolution equations 



^-^I-^-^ + f^ + ED + g 

S+ = H[-J(JV2 - iV 3 ) 2 - S + (3 + JL) + _|_(S 2 + 5| - 25?)] 

• „M-N% ro H. 4 7 r(S 2 -Sf) 1 
E_ = g 3 2 - 3 + 77? )E_ + — v 2 - 3; 
2v/3 # 2 \/3i7 2 

iV 2 = -iV 2 ff (-2E+ - 2V3E_ + 1 + — ) 
N 3 = -N 3 H (-2E+ + 2\/3E_ + 1 + — ) 

In analogy to Bianchi II these equations can be compared to (6.9)-(6.10) of [35] setting N\ to 
zero and using the definition of q (jT8]) . 



5 Special solutions 

In this section we present some special solutions which will be important, since we will show that 
the late time asymptotics of the Bianchi types considered behave like them in a sense which will be 
specified later. We start with the Kasner solution which is the general Bianchi I vacuum solution 
to motivate also the concept of generalized Kasner exponents. 
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5.1 The Kasner solution and generalized Kasner exponents 

The Kasner solution [T5] is the general Bianchi I vacuum solution, thus all components of the 
energy-momentum tensor and the scalar curvature R vanish. From the constraint equation one 
obtains: 

Y? + + E 2 . = 1 

which is known as the Kasner circle. The metric components are: 

3ii =diag(^,^,^) 

where pi, p 2 and p 3 satisfy: 

P1+P2+ P3 = 1 

Pi + p\ + pi = 1 

One can easily compute that the Hubble variable is H = 

For more general spacetimes let A^ be the eigenvalues of k%j with respect to <?jj, i.e., the solutions 

of: 

det(fc;- - A<5$) = (41) 

We define 

A, 

as the generalized Kasner exponents. They satisfy the first but in general not the second Kasner 
relation. 

5.2 The Collins-Stewart solution 

Another special solution which will play an important role is the Collins-Stewart solution ([5], p. 
430) with dust (7 = 1) which has Bianchi II symmetry: 

gcs = diag(2i, (2i) 3 / 2 , {2tf' 2 ) 

The Hubble parameter is H = ft" 1 and the energy density iirpcs = \t~ 2 ■ The values of the 
variables which have been introduced previously are: 

1 15 3 

S +=8 5 S -= 0;O =16 ;iVl = 4 

5.3 The Ellis-MacCallum solution 

In the case of Bianchi VIo there is a dust solution with diagonal metric discovered by Ellis and 
MacCallum (0, pp. 124-125): 

g E M = diag(i 2 , t 1 , t 1 ) 

arameter is H = |i 
8ttpem = t~ 2 is different. Here the values of the introduced variables are: 



The Hubble parameter is H = |t 1 as in the Collins-Stewart solution, but the energy density 



£+ = -I; E_ = 0; iVx = 0; N 2 = -N 3 = |; O = | 
The generalization of this solution to different values of 7 is called Collins solution. 
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6 Einstein-dust with small data 



6.1 Einstein-dust system 

Let us present briefly the Einstein-Euler system. We will consider the isentropic case where the 
matter fields are the energy density p and the four- velocity u a which is a unit timelike vector. The 
equation of state P = f(p) relates the pressure P with the energy density. The energy-momentum 
tensor is: 



T a = (p + P)u a u [3 + Pg a p 

and the equations of motion are equivalent to the condition that the energy-momentum tensor is 
divergence- free. The Einstein-dust system is then obtained via the condition P = 0. It can be seen 
as a very singular solution of the Einstein- Vlasov system. Formally the system can be obtained 
from the Einstein- Vlasov system choosing / to be of the form: 

f(t,x a ,p a ) = \u \\g (i) r i p(t,x a )S( P a - u a ) 

where uo is obtained via the mass shell relation. The relation of the Einstein-dust system to the 
Einstein- Vlasov system is in general subtle and we refer to [22] for more information on that. Here 
we will look at the special solutions of the corresponding Einstein-dust systems in order to obtain 
some intuition about the Einstein- Vlasov system. It is easy to see that the special solutions are 
equilibrium points. The stability of these equilibrium points has already been studied (see for 
instance [38 ). For the case of Bianchi I there even exist a general expression [(11-1.12) of [13]], 
where one can see that isotropization occurs. Actually for all the Bianchi cases we study here, 
Liapunov functions have been found, such that besides the stability also the global behaviour is 
known. In this section we will start dealing with estimates. C will denote an arbitrary constant 
and e a small and strictly positive constant. They both may appear several times in different 
equations or inequalities without being the same constant. 



6.2 Linearization of Einstein-dust around Collins-Stewart 



Let us look at the stability of the Collins-Stewart solution with dust. For the Collins-Stewart 
solution we have (S + = §,£_ = 0,iVi = §) which is an equilibrium point of the system (j38|) - (|40[) 
with S = 0. Let us translate the equilibrium point to the origin by introducing the variables 
t + = £+ - i £_ = £_ and Ni = JVi - §. The linearization is: 




93 





63 


64 


128 





3 
2 





87 







52 


— 64 




The variable £_ decouples and we obtain: 

£_ = £_(T )e-5( T - r «) 



The rest of the system: 



has eigenvalues 



S+\ _ Z-93/64 63/128\ (t+ 
N J ^-87/32 -3/64 J ^ 



Al/2 



(42) 
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Translated to our time variable means that the expected estimates for the Vlasov case are: 

E + -i = 0(r*) 

Whether this is true we do not know at this point. We will start proving estimates in the next 
chapter. Here we have obtained these estimates just in order to get a hint about the non-linear 
behaviour. For instance it could be the case that the variable £_ has the same decay as the other 
variables. However these estimate will turn out to be not sufficient. In addition to assume that we 
are close to the special solutions, we will have to assume that we are close to the dust case. This 
is done via a momentum bound. We have a number (different from zero) of particles at possibly 
different momenta and we define P as the supremum of the absolute value of these momenta at a 
given time t: 

P(t) = SU P {H = { g ab P aPb) L -\f{t,p) ? 0} 

A bound on that quantity can be used for estimates on S/H 2 as we show now. Consider an 
orthonormal frame and denote the components of the spatial part of the energy-momentum tensor 
in this frame by S a b- The components can be bounded by 

Sab < P 2 {t)p 

so we have that 

- < 3P 2 (43) 
P 

Since the Ricci scalar is non-positive it follows from © using the trace-free part of the second 
fundamental form that: 

167T/9 = 6H 2 + R- a ab a ab 

Thus we obtain: 

167T/9 < 6H 2 

Thus from [43] we obtain that: 



AttS 3 



< -P 2 . 



3H 2 ~ 2' 

For V due to ([25} and introducing the Collins-Stewart solution we have: 



V = -t-'g^V 2 - f^ig^V 2 + g 33 V 2 ) 
We see that V < —t~ l V which implies that the following holds: 

6.3 Linearization of Einstein-dust around Ellis-MacCallum 

For S = we have: 

S' + = -i(iV 2 -iV3) 2 -S + (2- g ) (44) 
N 2 - N 2 

E'_ = ^- ? p-(2-g)S_ (45) 

= N 2 (2E+ + 2 V3E_ + q) (46) 
= 7V 3 (2S+ - 2\/3£_ + q) (47) 
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with q=\ - 2i(iV 2 - N 3 ) 2 + §(5^ + E 2 „). 

For the Ellis-MacCallum solution we have (S + = — 4,£_ = 0,^2 = —N 3 = |) which is an 
equilibrium point of the system (|34" P -(|4"7 ]) . Introducing £+ = £+ + \, £_ = £_, iV 2 = N 2 - | 
and N 3 = N 3 



|. The linearization is: 



The eigenvalues are: 



PA 


( 


£_ 










\ 



21 
16 



15 
16 
_15 
16 



_3 

IV3 



_ 15 
32 

4 

_ J_ 
32 
_3_ 

32 



15 
32 

4 
J_ 

32 
__3_ 

32 





PA 




£_ 






) 


\N 3 J 



A1/2 = --(l±i\/3) 
A 3 /4 = -7(1 ±») 



(48) 



Translated to our time variable the expected estimates for the Vlasov case are: 

1 



N 2 
N 3 



4 
£_ 
3 

~ 4 
3 

+ 4 



= o(r 

0(t" 
0(t" 







With the same procedure as in the Bianchi II case, using now the Ellis-MacCallum solution we 
arrive at: 



p = o(H) 



7 The bootstrap argument 

The argument which will lead us to our main conclusions is a bootstrap argument, a kind of 
continuous induction argument. The argument will work as follows (see 10.3 of |28) for a detailed 
discussion). One has a solution of the evolution equations and assumes that the norm of that 
function depends continuously on the time variable. Assuming that one has small data initially at 
to, i.e. the norm of our function is small, one has to improve the decay rate of the norm such that 
the assumption that [to, T) with T < 00 is the maximal interval on which a solution with bounded 
norm corresponding to the prescribed initial data exists would lead to a contradiction. This is a 
way to obtain global existence for small data. In our case global existence is already clear but if 
the argument works we also obtain information about how the solution behaves asymptotically 
which is our goal. The interval we look at is [to,ti) and we will present the estimates assumed 
in the following for the different cases. All prefactors on the right hand side are positive and as 
small as we want. 



7.1 Bootstrap assumptions 

A first task is to find the suitable bootstrap assumptions. We choose a slightly slower decay for the 
anisotropy and the curvature variables than in the linearized cases with the hope that using the 
central equations, we are able to obtain the same decay as in the linearized case. For the estimate 
of P we start with a slower decay than the ones obtained in section [5] as well. The assumption 
of small data here is in the sense that our solutions are not "far away" from our special solutions. 
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In general to improve an estimate the corresponding evolution equation will be integrated. The 
assumptions made for the different Bianchi cases exclude the vacuum case, since the values of £1 
due to the constraint equation are near the corresponding values of O of the special solutions, thus 
far from being zero. 

7.1.1 Bootstrap assumptions for Bianchi II 

A+(l + t)~i 

A c (l + t)~i 

Amii + tyi 

7.1.2 Bootstrap assumptions for Bianchi VIq 





4' 


< 


A+(H 




3 
8 






< 


A_(H 


-ty 


3 
8 


\N 2 


4 1 


< 


A cl {l- 




3 
8 


\N 3 


+ 4 1 


< 


A c2 (l- 


M)- 


3 
8 




P 


< 


A m (l- 


vty 


1 

" 3 




p 



7.2 Estimate of the mean curvature 

The first variable we estimate is the trace of the second fundamental form or equivalently the 
Hubble variable. Let us rewrite (|2"U1) : 

d t (H- 1 ) = ^ + D (49) 

with 

Integrating (J^DJ) and since to = §-ff _1 (io) (this choice was made in section p.4[) ): 

H (t) = _L_ = V 1 



with 



Now for Bianchi II: 



and for Bianchi VIn: 



k + i 3 i + iit- 



I=[ D(s)ds 

Jtn 



2 K + ' 3H 2 24 
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D 



1 Ar , 2 

3^-24^-^ 



It turns out that in all cases D is small, in particular from the different bootstrap assumptions we 
obtain for for Bianchi II and VIq respectively: 



with 



\D\ < £2/3(1+0^ 



e 2 = C(A + +A 2 _+A c + A 2 m ) 

e 3 = C{A+ +A 2 _+ Aa. + A c2 + A 2 m ) 



We arrive at: 



'-r 1 ! = 0(e 2/3 t-i) 



The results for the Hubble variable is 



H=-r 1 (l + 0(e 2/3 rl)) 



(50) 



(51) 



We also obtain an estimate for the determinant using the estimate of H and integrating © in 
both directions. 



C(t )t 4 - e < g(t) < C(t )t 



4+e 



(52) 



7.3 Estimate of the metric 

Consider the following equation in the sense of components: 

g ab = th ab 

In particular we will consider the components g 22 and <? 33 , which means that for Bianchi II 2 = | 
and for Bianchi VIn - = 1. We will show that: 

u n 



±(t^g ab ) = r^- 1 3 a6 (-7 + - + ^- b t)< -nr^-n ab 

at q g a ° 



^9" 



with rj positive with the help of the bootstrap assumptions and choosing 7 in a suitable way. This 
means then that: 



|rV)<o 



which implies what we wanted to show: 



g ab (t) < t^ g ab (t )t-^'. 

p 

For the covariant components one can do the same by defining g a b = t q g a b- One obtains: 



(53) 
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at q g a b 

For the last step one can actually use the same 7 as for the contravariant components since 
9abg ab = —gab9 ab - In other words once is shown, we also have: 



9ab{t) < t 9 9ab{ta)t" 7 
From the definitions made one can obtain: 

g 11 =2. 9 n tf(-l + 2I] + ) (54) 



2g 22 H{-\ - £+ - \/3S_) (55) 
g „„ = 2.g 33 J ff(-l-S+ + V3S_) (56) 



3 22 

^33 



Then we have with (|5lT) for the components g 22 and g 33 : 



V 



i] = j + 2Ht(l + E + ±V3E-) - 

= 7+|(l + 0(e 2/3 t-i))(l + S + ±V3S_)-^ 
6 q 

In both Bianchi II and VLj: 

|(1 + S+) -5 = 0(^(1 + *)-*) 

which enables us to choose 7 in such a way that 77 is positive. Different values of S + correspond to 
different exponents in the components of the metric. Using the estimates of g 22 and g 33 we obtain 
then the estimate for the other component of the metric <?n via the estimate of the determinant. 
We could also proceed directly from ([5"4"|) . 

Summarizing this means that asymptotically up to a positive constant which depends only 
on to the components (and their inverses) of the metrics gn for Bianchi II and gvi f° r Bianchi 
VIo have the same decay up to an e as the corresponding components of the Collins-Stewart and 
Ellis-MacCallum solution respectively: 

C(* )*- e < — < C(t )t+* 
9cs 

C(t )t- e < ^ < C(t )t+< 
gEM 

7.4 Estimate of P 

We can express the derivative of the metric as follows: 

g»f = 2H{Y b a - S b a )g a f 
It follows from and using (j2"6]) : 

V = g bf V b V f = 2ff(E» ~ 5 b a )g a fV b V f = 2H{^\g 11 V 2 + ^\g 22 V 2 + Z 3 g 33 V 2 ) - 2HV 
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The maximum of Ej, E| and E§ is for Bianchi II and VIo equal to \ + 0(t « ). Thus: 

V < 2HV(— + et-*) 
Using now the estimate of H and integrating : 

v<v(t )(t/t )- 1+t 

from which follows: 



i— e 

Choosing P(to) < A m t§ we arrive at: 

P < A m t-i +e 

which is an improvement of the bootstrap assumption and which has the consequence that: 



— < Ct~ 1+t 
H 2 ~ L 



(57) 



7.5 Closing Bianchi II 

Until now we have estimates for H and for P in the interval [to, ti). We need to improve the other 
variables. Although Bianchi II and Bianchi VIo are more complicated, the main argument will be 
the same as in Bianchi I |21) . 

7.5.1 Estimate for E_ 

Case I If |E_| < A_(l + t)~ 1+e holds there nothing more to do, since this is a better estimate 
then the one assumed. 

Case Ila Assume now E_ > A_(l + t)~ 1+t . Define t<x as the smallest number not smaller than 
to with the property E_ > A_(l + t)~ 1+e . Since we are assuming that E_ > we can divide (|36l) 
by E_: 

§- = ^[-(S+^+El 1 ^^-^)] 

With l|5TI] ) . ((51 11 . tlie fact that s 2 - s 3 < S, and our assumption: 

|^ < -^(l-O (58) 

where £ is as small as we want. This variable £ contains a term of type AZ 1 A 2 n , but A_ and A m 
can be chosen independently as small as needed. Also the e coming from S has to be chosen bigger 
than the e coming from E_. Note that it becomes clear here why we had to improve P to arrive 
at (|57p . Integrating ((55)) between ti and f 2 : 

E_(tO < E_(t 2 )t^ 1+ « 

Assume i 2 = then: 

E_(to < s_(to)t;- c tr 1+e < A_t- l+i 

since E_(to) can be chosen in such a way that the last inequality holds. 
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Case lib If £2 > to then by continuity £_(ia) < A-(l + 1) 1+£ which means that: 

E_(ti) < + t) £ -^ 1+5 < 

if e is chosen to be smaller than £. The argument for the case that E_ is negative is the same, 
just define E_ = — E_ and use if?f — Sf < 5. This means that we could improve our bootstrap 
assumption to: 

|S_(ti)| < 

7.5.2 Bootstrap assumptions for the other time variable 

We have found that for the estimates in the following section it was useful, although not essential, 
to use the other time variable r. Using the estimate of the Hubble variable (j5Tj) in the definition 
of t (flu)) we have: 

T-T = j \t- X {l + 0{drl))dt 

After integrating and observing that tq and to are constants we arrive at: 



-T + To + C 

where £ is small: £ = 0(e(t~« + t 8 )). So the bootstrap assumptions can be translated to the 
time variable r. We obtain: 



E+| 


< 


CA+e 




E-| 


< 


CA_e 




|JVk| 


< 


CA c e~ 


A- 


P 


< 


CA m e 





Since we have an estimate of H in both directions one can go also back from an estimate in terms 
of r to an estimate of t just by a multiplication by a constant which will not be relevant. 

7.5.3 Estimate for £+ and Ni 

Define 

Nj 11 \N ly 

where Mn is the matrix of eigenvectors of the linearized system (|^2"]l . Then we have: 

l)'H(_)jf)(l:H^>C 

since 0(£+ + N? + E 2 . + P 2 ) = O^^e- 1 "). Multiplying the first equation by E + and the second 
by N± and adding both we obtain: 
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Let us assume now that: 

t? + + N* > (Al+A 2 c )e(-i+V T 

This implies: 

Z 2 + + Nl > C(A\ + Al)e { -i + ^ 

From which follows that: 

t\+Nl < (E^(r )+iV 1 2 (r )) e ^i +e )( T - T «) 
< C(^(t )+iV 1 2 (to))(f)- 1+t 

TO 

or: 

t\ + N? < C(£ 2 + (t Q ) + JV?( fo ))(±)-i+« 

to 

Making now the same argument as in the end of the estimate of £_ we arrive at improved estimates 
for £+ and N r : 

|£+| < t + (t Q )t-^ 

1^1 < iV 1 (t )t-i+ £ 

i.e. for S + and N\: 

< A + (i + i )-i+ e 
|iVi-||< A c (i+t)-=+« 

We have closed now the bootstrap argument. Note that for this last improvement of the estimates 
S + and N% we did not use the improved estimates for £_ and P. 

7.6 Closing Bianchi VI 

This case is analogous to Bianchi II. The bootstrap assumptions with the variable t read: 





< 


CA+e- 


16 T 




< 


CA-e~ 




n 2 \ 


< 


CA cl e~ 




n 3 \ 


< 


CA c2 e~ 




p 


< 


GA m e~ 





In this case what remains are the estimates for £+, N 2 and N3. In terms of the transformed 
linearization 









iV 2 


= M VW 




\N 3 J 




\N 3 J 
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where Mvio is the matrix of eigenvectors of the linearized system (J2HJ) we have: 



PA 




/ 1 


-V3 


°\ 






3 




1 














1 


-1 


^2 


\N 3 J 






1 


1/ 





(A 
1 

1 

vv 



^2 , M2 + ^2 + £2 

same procedure at 



since 0(Y? + + N$ + N$ + Y, 2 _ + P 2 ) = 0(A 2 n e~ T ). As in the Bianchi II case, we arrive with the 



(fl + + + £ 2 _)- 1 0(A 2 m e- T ) 



|[log(iV| + iV 3 2 )] = -| + (# a + iV 3 )(jvf + N 2 r l O{A 2 m e-T) 



and this means that: 



±[lo g (± 2 + + t 2 _)]<-l + ee^-^ 



and a similar expression for N 2 and ./V3 such that in the end we arrive at the estimates we wanted 
to obtain. In this case as well it was not necessary to use the improved estimate of P. 



7.7 Results of the bootstrap argument 

Since we have improved all estimates we have closed the bootstrap argument. Let us summarize 
the results obtained in this chapter in the following propositions: 

Proposition 1. Consider any C°° solution of the Einstein- Vlasov system with reflection and 
Bianchi II symmetry and with C°° initial data. Assume that \T, + (to) — g|> |S_(to)|, \Ni(to) — || 
and P(to) are sufficiently small. Then at late times the following estimates hold: 

V 1 (l + 0(t"5+ £ )) 



H(t) 
1 

~ 8 
£_ 
3 

~ 4 
P(t) 



3 

o(t- 1+£ 
o(t-^ +e 
o(t-? +e 



Proposition 2. Consider any C°° solution of the Einstein- Vlasov system with reflection Bianchi 
VIq symmetry and with C ,0 ° initial data. Assume that |£+(to) + ^|, |£-(^o)|; I -^2 (^o) — §|> I -^3(^0) + 
1 1 and P(to) are sufficiently small. Then at late times the following estimates hold: 



H(t) 


= ^(l+orh) 




= 0(H+ £ ) 




= o(t-^ +e ) 




= o(t-^ +e ) 


N a + \ 


= o{r^ +t ) 


P(t) 


= o(t-i +e ) 



In the next chapter we will improve the estimates such that we can get rid of the e. However 
the results stated here in this section represent in fact the core of our results. 
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8 Main results 



8.1 Arzela Ascoli 

Until now we have obtained estimates which show that the decay rates of the different variables 
are up to an e the decay rates one obtains from the linearization. We want to use the Arzela- 
Ascoli theorem. We will show the boundedness of the relevant variables and their derivatives. 
The variables E+, N±, N 2 and N 3 corresponding to the different Bianchi cases are bounded 
uniformly due to the constraint equation. In particular: 

Y? + + S 2 . < 1 
N 2 < 12 
(N 2 - N 3 ) 2 < 12 

Note that iV3 is negative. The Hubble variable H and its derivative is bounded as one can see 
from the estimates and [30] We have also obtained with the bootstrap argument that P, which 
is non-negative, decays which means that S/ H 2 is bounded. From the estimates obtained it is 
clear that g ab and its derivative are bounded. Now having a look at the central equations we see 
that the derivatives of E + , N±, N2 and N 3 are also bounded uniformly. If we can bound the 
derivative of S also the second derivatives of E+, N±, N 2 , N 3 and H are bounded. For this 
purpose it is convenient to express the components of the energy momentum tensor in terms of 
integrals of the covariant momenta: 

S ab g ab = J f(t, P ) PaPb g ab (l+g cd p cPd )-ig- 1 2 dpid p2d P3 
= J f(t,p)V(l + V)-ig-U Pl dp 2 dp a 

The only term of the time derivative of S which could cause problems is the time derivative of the 
distribution function, since V and g ab can be bounded by V and g ab respectively and we know that 
S itself is bounded since S/H 2 is. The term with the time derivative of the distribution function 
can be handled with the Vlasov equation: 

J f(t,p)V(l + V)-?g-id Pl dp 2 dp 3 
= - J (pVC*P 6 P<i|£v(1 + V)-?g-?dp 1 dp 2 dpz 



Integrating by parts we obtain a term which can be bounded by S. Note that the momenta 
grow in the worst case with t 1 and that po is also bounded from below since the particles are 
assumed to have mass. Now all the relevant quantities are bounded. Let {t n } be a sequence 
tending to infinity and let (E_)„(t) = E_(i + t n ), (S+)„(t) = E+(t + t B ), (JVi) n (t) = iVi(f + i„), 
(N 2 )n(t) = N 2 (t + tn), (N 3 ) n (t) = N 3 (t + tn), H n {t) = H{t + t n ) and S n (t) = S{t + t n ). Using 
the bounds already listed, the Arzela- Ascoli theorem can be applied. This implies that, after 
passing to a subsequence, (E_)„, (E+) n (Ni) n , (N 2 ) n , (N 3 ) n , H n and S n converge uniformly on 
compact sets to a limit (E_)oo, (E+)oo (iVx)oo) (-^2)00, (N 3 )oo, Hoo and 5*00 respectively. The first 
derivative of these variables converges to the corresponding derivative of the limits since we have 
been able to bound the derivative of S in the last section. Going to this limit it is easy to see that 
the variable D of equation is zero and consequently: 

Una = — t ^ 

3 

From (|M)) - ([55)) we see that for Bianchi II and VIo we obtain the optimal decay rates for the metric 
and for its derivative. This implies that we obtain the optimal decay rates for P. Since S/H 2 is 
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zero asymptotically we obtain the same estimates for £_, E+, Ni, N-2 and 7V3 as in the Einstein- 
dust case. Introducing this estimates in (f4U)h we also obtain the optimal estimate for H. Let us 
summarize the estimates. 

8.2 Optimal estimates 

Theorem 1. Consider any C°° solution of the Einstein- Vlasov system with reflection and Bianchi 
II symmetry and with C°° initial data. Assume that |£+(fo) — g|, |£_(to)|, \Ni(to) — || and P(to) 
are sufficiently small. Then at late times the following estimates hold: 

H(t) = H t -i(i + (t-*)) 
£_ = O^" 1 ) 

JVi-| = O(i-i) 
P(t) = 0(t"*) 

Theorem 2. Consider any C°° solution of the Einstein- Vlasov system with reflection Bianchi VIq 
symmetry and with C°° initial data. Assume that |£+(fo) + \\, |£_(io)|? I -^2 (^o ) — f| ; l-^3(^o) + f| 
and P(to) are sufficiently small. Then at late times the following estimates hold: 

H(t) = J-\l + 0(t-i)) 

E + + i = 0(H) 
£_ = 0(i"5) 
N 2 -\ = 0(t~i) 

n 3 +'± = O(t-i) 
P(t) = 0(t~i) 

For the cases Bianchi II and VIo we are also able to obtain the optimal estimate for the metrics: 

Corollary 1. Consider the same assumptions as in the previous theorem concerning Bianchi II 
and VIq respectively. Then 

gu = tdiag{K 1 ,t 1 ' 2 K 2 ,t 1 ' 2 K 3 ) 
gvi = tdi&g(tK 4 , K 5 , K e ) 

with K n = C n + 0(t~i) and where C\-C§ are independent of time. The corresponding result for 
the inverse metric also holds. 

We see that the error in the metrics comes from the error in £ + . 

8.3 Kasner exponents 

From (jT3J) we see that the eigenvalues (|4ip of the second fundamental form with respect to the 
induced metric are also the solutions of: 

det(aj - [X-\k]Sj) =0 
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Let us define the eigenvalues of cr^ with respect to gij by A^, we have that: 

A 4 = Aj - -k 

Note that Ei(Aj) 2 = a a b<r ab . In the cases Bianchi II and VIo since everything is diagonal the 
Kasner exponents are easy to calculate. Using the optimal estimates for £ + , £_ and H and the 
fact that the sum of the generalized Kasner exponents is equal to one, we finally arrive at the 
generalized Kasner exponents for Bianchi II which are (j, §, §) and for Bianchi VIo (3, \ , \ ) in 
both cases up to an error of order 0(t~^). Let us summarize these result: 

Corollary 2. Consider the same assumptions as in the previous theorem concerning Bianchi II 
and VIq respectively. Then: 

Pn = Pes + 0(t~i) 
Pvio = Pem + 0{t~i) 

We see that the error in the Kasner exponents comes from the error in H . 
8.4 Estimates of the energy momentum tensor 

Before coming to the estimates of the energy momentum tensor we show that for Bianchi II V2 and 
V3 become constants, something similar can be done for Bianchi VIo- Define E = g 22 V 2 2 + g 33 V 2 , 
then: 

E = g 22 V£ + g 33 Vf < 2H(-1 - S+ + V3\Z_\)E 

Integrating 

log[E/E(t )} = -^logt/to + 0(e(t-? 
We have the following inequality for E: 

E < Ct~i 

Since the components of the metric g 22 and g 33 tend to the corresponding components of the 
Collins-Stewart solution we see that V2 and V3 become constant asymptotically The same is true 
in the case of the Ellis-MacCallum solution. Now since f(to,p) has compact support on p, we 
obtain that there exists a constant C such that: 

/(*,*») = \Pi\>C 

Let us denote by p the momenta in an orthonormal frame. Since f(t,p) is constant along the 
characteristics we have: 

|/(i,p)|<||/o||=sup{|/(i ,p)|} 
Putting these facts together we arrive at the estimates which we summarize in the following: 

Corollary 3. Consider the same assumptions as in the previous theorem concerning Bianchi II. 
Then 

P = Pcs(l + 0(t-i)) 
S- tJ < C\f \t- 3 
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Corollary 4. Consider the same assumptions as in the previous theorem concerning Bianchi VIq 
respectively. Then 

P = pEM(l + 0(t~i)) 
Sij < C\f Q \t~ 3 

The error in the energy density comes from the error in H . 

Remark From the corollaries one can estimate the quotient Sij/p which is 0(t~ l ). That 
this quotient vanishes asymptotically means that the matter behaves as dust asymptotically as 
expected. 

9 Conclusions and Outlook 

The results concerning Bianchi II generalize the results obtained in [3D]- For Bianchi VIo even 
for the reflection symmetric case there is no analogous previous result. The reason is that it is 
not compatible with the LRS-symmetry. Thus our result concerning Bianchi VIo shows clearly 
that the methods developed are powerful in the sense that one can obtain results which where 
out of reach with the techniques developed until now. However in contrast to |30| our results are 
restricted to the case of small data. In order to remove this assumption the Liapunov functions 
discovered for the fluid model could be helpful. In the case of Bianchi II the future asymptotics 
are known globally even in the tilted case [15] . 

We also hope to extend these results to the case without reflection symmetry in a future 
publication. 

In our argument we have used the Arzela-Ascoli theorem, but only at the end. Thus there 
exist a lot of estimates where one has control over the constants involved. Probably this could 
help for a numerical analysis of the Einstein- Vlasov equation which is quite difficult. 

Another path of generalizing our results could be the extension to higher dimensions. For the 
vacuum case geodesic completeness was shown for some homogeneous models in higher dimensions 

. The work on homogeneous spacetimes in higher dimensions may also shed some light on the 
inhomogeneous case in four spacetime dimensions. 
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